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We apply the Nose-Hoover thermostat and three variations of it, which control different combina- 
tions of velocity moments, to the periodic Lorentz gas. Switching on an external electric field leads 
to nonequilibrium steady states for the four models with a constant average kinetic energy of the 
moving particle. We study the probability density, the conductivity and the attractor in nonequi- 
librium and compare the results to the Gaussian thermostated Lorentz gas and to the Lorentz gas 
as thermostated by deterministic scattering. 



I. INTRODUCTION 



In a system of particles under an external force a nonequilibrium steady state can be obtained by applying a 
thermostat [|]-f§. Deterministic and time reversible bulk thermostating is based on introducing a momentum de- 
pendent friction coefficient in the equations of motion. One type of this mechanism is the Nose-Hoover thermostat 
. It creates a canonical ensemble in equilibrium and yields a stationary nonequilibrium distribution of velocities 
in nonequilibrium. Another version, the Gaussian isokinetic thermostat |6|-|8|, leads to a microcanonical density for 
the velocity components in equilibrium and to a constant kinetic energy in nonequilibrium. Though the microscopic 
(N ■ dynamics of these thermostated systems is time reversible the macroscopic dynamics is irreversible in nonequilibrium 
M . This is related to a contraction onto a fractal attractor |lC|-|l2]] . 

Characteristic features of thermostated many particle systems, like a nonequilibrium steady state and a fractal 
attractor, have been recovered for a specific one particle system, the Gaussian thermostated Lorentz gas Jl0|,[l3] pS) . 
Ch ■ The periodic Lorentz gas consists of a particle that moves through a triangular lattice of hard disks and is elastically 
^T >' reflected at a collision with a diskQ It serves as a standard model in the field of chaos and transport, see e. g., P,[ui|. 
In contrast to many particle systems a one particle system reflects more strongly the properties of a thermostat. 
For the Gaussian thermostated Lorentz gas a complicated dependence of the attractor on the field strength results 

©mi- 

A second simple one particle system which has been much investigated concerning its chaotic properties is the Nose- 
Hoover thermostated harmonic oscillator. In contrast to the Lorentz gas, the dynamics of this system is generically 
noncrgodic 

Hll. However, one can obtain an ergodic dynamics for this system by the additional control of the 
| square of the kinetic energy Jl9[-^l[|. 

Recently an alternative thermostating mechanism, called thermostating by deterministic scattering, has been intro- 
duced for the periodic Lorentz gas ^^3|Q. This deterministic and time reversible mechanism is based on including 
energy transfer between moving particle and disk scatterer at a collision, instead of using a momentum dependent 
friction coefficient. It leads to a canonical probability density for the particle in equilibrium, and in nonequilibrium it 
keeps the energy of the particle on average constant. In Refs. 22 2jJ this model has been compared to the Gaussian 
thermostated Lorentz gas: In nonequilibrium one finds an attractor for this model which is similar to the fractal 
attractor of the Gaussian thermostated Lorentz gas, but in contrast to the Gaussian case the attractor is phase space 
filling even for high field strengths. For both models the conductivity is a nonlinear decreasing function with increasing 
field strength on a coarse scale. 

Based on its construction the method of thermostating by deterministic scattering is in fact closer to the Nose- 
Hoover thermostat. This motivates us to apply for the first time the Nose-Hoover thermostat to the periodic Lorentz 



*Note: Most of the figures are in poor quality output. The originals are many MB large. They can be obtained upon request. 
: A model almost identical to the driven periodic Lorentz gas, except for some geometric restrictions, is the Galton board, 
which has been invented in 1873 to study probability distributions Jl 7| ] . 
2 This mechanism has been applied later to a system of hard disk under a temperature gradient and shear E3. 
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gas. In Section II we introduce the Nose-Hoover thermostat and discuss some variations of it. In Section III we define 
the periodic Lorentz gas and the thermostats we will study. We investigate these models in equilibrium in section IV 
and in nonequilibrium in Section V. Conclusions are drawn in Section VI. 



II. THE NOSE-HOOVER THERMOSTAT AND SOME VARIATIONS 



In the following sections we consider a one particle system in two dimensions with the position coordinates q = 
(q x ,q y ) and the momentum coordinates p — (p x ,p y ). The mass of the particle has been set equal to unity. The 
equations of motion for the Nose-Hoover thermostat are then given by H 



q=p 

P = £- (P 

V 2 I 
•2T t 



The thermostat variable C couples the particle dynamics to a reservoir. It controls the kinetic energy of the particle 
p 2 /2 such that < p 2 >= 2T. This holds even in nonequilibrium as induced by an electric field e. r is the response time 
of the thermostat. Performing the limit r — > in Eqs. ([!]) approximates the Gaussian thermostat with £ = (e- p)/p 2 . 
In the limit t — * oo the friction coefficient approaches a constant, C = Co and the equations of motion are not time 
reversible anymore. The dynamics of this dissipative limit has been investigated in [25[ . 

A generalization of the Nose-Hoover thermostat to control higher even moments of p has been introduced by Hoover 
p6| . The moments are fixed according to the momentum relations for the Gaussian distribution. By such a more 
detailed control of the nonequilibrium steady state statistical dynamical properties, like ergodicity, can be improved, 
as has been mentioned for the Nose-Hoover thermostated harmonic oscillator in the introduction. In principle the 
method of the control of the even moments can straightforwardly be extended to the control of the odd moments. 
However, such a thermostat would involve an additional parameter to control the respective current of the subsystem. 
Thus, the corresponding reservoir would be more than a single thermal reservoir, which is physically not desirable, 
apart from the fact that the respective equations of motion would not be time reversible an ym ore. 

We briefly note that there exist other formal generalizations (^tJ, or modifications pj|p9| , of the Nose-Hoover 
thermostat in the literature. They have been critically reviewed in Rcfs. p(i|,pil|,f26[ . 



III. VARIATIONS OF THE NOSE-HOOVER THERMOSTAT FOR THE PERIODIC LORENTZ GAS 

The basic themostating method we investigate in this paper is the Nose-Hoover thermostat, Eqs. ([!]). In the follow- 
ing we introduce three variations of it which are time reversible and result in a dissipative dynamics in nonequilibrium. 
The first variation is the Nose-Hoover thermostat with a field dependent coupling to the reservoir, 

q=p 

Px = £x - (1 + £x)(Px 
Py = £y - (1 + £y)(Py 

which is obtained by including the factors 1 + e x , resp. 1 + e y , in Eqs. (0). Alternatively, these equations can be 
written by defining two field dependent friction coefficients, £ x = (1 + e x )( and £ y — (1 4- £ y )C, which are governed 
by £ x — (p 2 /2T — 1)(1 + s x )/t 2 and £ y = (p 2 /2T — 1)(1 + e y )/r 2 , respectively. It then becomes clear that for each 
momentum component there exists a separate response time of the reservoir which reads + e x , resp. r/yT + e y . 

The basic advantage of this thermostat is that the response times are now adjusted to the corresponding component 
of the field strength such that with increasing field strength the response time decreases. The standard Nose-Hoover 
thermostat Eqs. ([!]) is contained as a special case in equilibrium . 

The second variation goes back to Hoover [^6|. It includes a control of < p 4 >= 8T 2 , 

q=p 

p 2 

p = e- Cip- C, 2 —p 
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p 2 p 2 1 



C 2 -^(^-2)^. (3) 

As mentioned in the previous sections this variation can improve statistical dynamical properties, like ergodicity. 

The third variation controls p 2 and p 2 separately. This is performed by defining two independent reservoirs for the 
x- and y-direction, 

q=p 

Px — &x C,xPx 

Py = £ y ~ CyPy 
V 2 1 

c = (|-i)^ 



4 = (f-i)^- (4) 

v 

This variation more deeply intervenes in the microscopic dynamics by forcing the single components p 2 and p 2 
separately towards canonical distributions. However, in contrast to the previous variations a curiosity is hidden in it: 
At a collision the thermostat £ x , £ y is uncorrelated to the thermostated variables p x and p y , because p x and p y change 
at a collision whereas C, X ,Q V remain the same. Therefore the thermostat does not work efficiently. But we have not 
found any reflection of this curiosity in the macroscopic behavior. 

We study the dynamics of these models in one Lorentz gas cell with periodic boundaries, see Fig. 0(a). As the 
radius of the disk we take r = 1. For the spacing between two neighboring disks we choose w ~ 0.2361, corresponding 
to a density equal to 4/5 of the maximum packing density of the scatterers. A collisionless free flight of the particle 
is avoided for this parameter p0[ . The relevant variables of the dynamical system are defined in Fig. 0(b): (3 is the 
angular coordinate of the point at which the particle elastically collides with the disk, and 7 is the angle of incidence 
at this point. 

The equations of motion are integrated by a fourth order Runge Kutta algorithm with a step size of dt = 0.005 
between two collisions. The collision of the particle with the disk has been determined with a precision of 10 -7 . 
Unless declared otherwise the temperature is chosen to T = 0.5. 

IV. EQUILIBRIUM 

Inserting the initial condition (pojCo) = (V2T, 0) in Eqs. (0) with e = the velocity of the particle becomes a 
constant, p — \/2T, that means the Nose-Hoover thermostat does not act in the Lorentz gas and the dynamics is 
microcanonical. For other initial conditions one observes in computer simulations that p 2 and £ oscillate periodically 
and that the dynamics of this one particle system is nonergodic. 

A stability analysis confirms the numerical results: The Nose-Hoover thermostated equations of motion Eqs. (0) 
can be reduced for e — to 



p 2 = -2(p 
V 2T 



C=(^-l)-. (5) 



Eqs. (||) are also valid at the moment of a collision, because p 2 and £ are not changed by a collision. The fixed point 
of Eqs. (||) is (p 2 , C) = (2T, 0) with the eigenvalues A1/2 = ±a/— 4T/t 2 and is thus elliptic. 

The additional control of < p 4 > destroys the microcanonical probability density but it is not sufficient to obtain 
an "exact" dynamics^ in the periodic Lorentz gas. Different initial conditions still lead to a different shape of the 
probability density g(p x ), as is shown in Fig. ||(a). 

In contrast the separate control of p 2 and p 2 leads to an exact dynamics in equilibrium corresponding to the 
canonical probability density g(p x ), as shown in Fig. 0(b). 



3 The dynamics is "exact" if every initial density of nonzero measure converges to the same stationary density 
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V. NONEQUILIBRIUM 

We now apply an external electric field e parallel to the x-axis. The Nose-Hoover thermostat and the related 
models then lead to well defined nonequilibrium steady states with constant average energy of the particle. 

A. Probability density g(p x ) 

The probability density g(p x ) for the Nose-Hoover thermostat for e = 0.5 is presented in Fig. ||(a). For r 2 = 0.01 
the density shows some remains of the deformed microcanonical density of the Gaussian thermostat, whereas for 
t 2 = 1 and t 2 = 1000 the density becomes similar to the density of thermostating by deterministic scattering, which 
is related to a canonical distribution f p2p3j ] . 

Fig. ||(b) shows g(p x ) for the three variations of the Nose-Hoover thermostat. We have chosen here T = 0.60029 
which corresponds to the temperature in the bulk for thermostating by deterministic scattering at a parametric 
temperature of T = 0.5. The density of the Nose-Hoover thermostat with field dependent coupling to the reservoir 
is very close to the density of thermostating by deterministic scattering. The density of the Nose-Hoover thermostat 
with separate control of p\ and p y looks like a superposition of the densities of the Nose-Hoover thermostat for small 
and for large r. 

In all models the mean value of g(p x ) is positive, indicating a current parallel to the field direction. 

B. Conductivity 

The conductivity a =< p x > je for the Nose-Hoover thermostat is shown in Fig. ||. For r 2 = 0.01 the curve is 
very similar to the conductivity of the Gaussian thermostated Lorentz gas fy| . For r 2 = 1000 the curve is more 
stretched along the e-axis and globally not decreasing anymore. In contrast the conductivity as obtained from 
thermostating by deterministic scattering J2^,^3| is a globally decreasing function. According to the Einstein relation, 
in the limit e — > a should approach the equilibrium diffusion coefficient D of the periodic Lorentz gas, which for 
w = 0.2361 has the value D w 0.21 [[IB). This is hard to see for r 2 = 1000, because for e — > the probability density 
changes drastically from a smooth, canonical like density to a non-smooth density. It is also difficult to see any linear 
response in computer simulations, as has already been discussed for the Gaussian thermostated Lorentz gas and for 
thermostating by deterministic scattering in Ref. p3] . 

C. Attractor 

Fig. H shows the Poincare section of [fi, sin(7))R at the moment of the collision for the Nose-Hoover thermostat, 
for the three variations and for thermostating by deterministic scattering. Again, we have chosen T = 0.60029 to 
compare the results with thermostating by deterministic scattering. For all models the structure of the attractor 
is qualitatively the same as the structure of the fractal attractor obtained for the Gaussian thermostated Lorentz 
gas J[(j| ]. However, the fine structure varies with the models and with the response time r. For the Nose-Hoover 
thermostat with t 2 = 0.01, see Fig ||(a), the structure is most pronounced, whereas for the control of p\ and p y 
separately, see Fig. 0(e), the structure is least visible. 

D. Bifurcation diagram 

The angle (3 at the moment of the collision is presented as a function of the field strength for the Nose-Hoover 
thermostat in Fig. |(| For all three values of r 2 the attractor is phase space filling for small field strengths e < 1.3 
and contracts onto a periodic orbit with increasing field strength. For r 2 = 0.01 the scenario is similar to the one 
of the Gaussian thermostated Lorentz gas Jlj|0. For r 2 = 1 the scenario looses its richness, but it gets a bit more 



4 Because the symbols for the angles vary in the literature we mention again that the angle j3 gives the location of the collision 
relative to the field direction and 7 is the angle of incidence. 
5 In Fig. 10 of [Q the angle of flight after a collision is plotted and the field is parallel to the negative x-axis. 
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complicated again for r 2 = 1000. In contrast to the Nose-Hoover thermostat, the attractor of thermostating by 
deterministic scattering remains phase space filling even for large el 2^,^3| . 

The bifurcation diagram in the dissipative limit r — » oo of Eqs. (|l) with a constant friction coefficient ( c shows an 
inverse scenario to Fig. ||, as is presented in Fig. |j. For small e the trajectory is a so-called creeping orbit ^5|, then 
it changes to a periodic orbit, and for large e the attractor gets phase space filling. By increasing £ c the strength 
of the dissipation increases with the consequence that the onset of chaotic behavior starts at higher field strengths. 
With respect to the numerics we remark that related to the different shapes of the attractors by varying £ c in the 
dissipative limit, the duration of the transient behavior of the Nose-Hoover thermostat grows drastically for r — > oo. 

Figs. p[]l0| show the bifurcation diagrams for the three variations of the Nose-Hoover thermostat. In general one 
observes that even for high field strengths chaotic regions appear, in contrast to the Nose-Hoover thermostat. 

For the variation with the additional control of p A the attractor covers a bounded /3-interval for e > 3. For these 
field strengths the trajectory is a creeping orbit. 

Fig. |l^(b) depicts the attractor for the variation with separate control of p 2 and p 2 under different response times 
for the x- and y-direction. Since the field acts in the x-direction we have chosen a strong coupling, t 2 = 0.1, for 
the x-direction and a weak coupling, t 2 = 1000, for the y-direction. Up to a field strength e m 6.5 no periodic 
window has been found. The bifurcation diagram for these parameters most strongly deviates from the Nose-Hoover 
and Gaussian thermostated Lorentz gas and is closest to the bifurcation diagram of thermostating by deterministic 
scattering. However, in contrast to thermostating by deterministic scattering the attractor is more concentrated 
around (3 ~ 7T. 

We detected a numerical problem for the second and the third variation at several values of r. One observes that 
sporadically after large time intervals there appears a creeping orbit with a very low velocity of the particle, which 
is difficult to handle numerically. Whether this creeping orbit is the stationary state could probably be clarified by 
calculating Lyapunov exponents flq] . 



E. Thermodynamic entropy production and phase space volume contraction 

A characteristic property of the Nose-Hoover thermostat as well as of the Gaussian isokinetic thermostat is that 
the thermodynamic entropy production is equal to the phase space volume contraction rate |J. This equality can as 
well easily be verified for the variation with the additional control of p A and for the variation with separate control of 
pi and p 2 y . 

On the other hand it does not hold for the Nose-Hoover thermostat with field dependent coupling to the reservoir. 
From Eqs. (||) one gets for the phase space contraction rate of this model — < divr >= (2 + e) < £ > where 
r = (q,p, C)- The precise relation between thermodynamic entropy production Std — £ < Px > /T and — < divr > 
for this variation is obtained by calculating the energy balance between subsystem and reservoir: 

2 

E = V — + Tt 2 C (6) 
is the total energy which in a noncquilibrium steady state should on average be zero, 

<f>- 



Inserting Eqs. (g) with e y = in Eqs. (0) leads to 

£x < Px > £x < plC > 



T T 



2 < C > • (8) 



Numerical simulations have shown that p 2 and £ are no independent quantities in noncquilibrium. If p 2 and C would 
be independent and equipartitioning would be fulfilled, i. e., < p 2 >= T, which is only the case in equilibrium, then 
Eq. (@) would lead to an identity between thermodynamic entropy production and phase space volume contraction. 
The results for — < divr > and for Std as obtained from computer simulations for this system are presented in Table 
I at different r and s. More details of the entropy production in this variation and in a Gaussian thermostat with the 
same property are discussed in Ref. |Rl| . 
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VI. CONCLUSIONS 



We have investigated the Nose-Hoover thermostat and three variations of it for the periodic Lorentz gas. All models 
are time reversible and lead to well defined nonequilibrium steady states with a constant average kinetic energy of 
the moving particle. 

As a typical characteristic of deterministic and time reversible thermostating mechanisms it has been confirmed 
that in nonequilibrium all these systems contract onto attractors similar to the fractal attractor of the Gaussian 
thcrmostated Lorentz gas. 

In equilibrium only the variation of the Nose-Hoover thermostat with separate control of p 2 x and p y leads to an exact 
dynamics with g(p x ) being canonical just like the corresponding density of thermostating by deterministic scatterin 
g- 

In nonequilibrium the attractor of the Nose-Hoover thermostat contracts onto a periodic orbit for higher field 
strength, analogous to the Gaussian thermostated Lorentz gas. However, the detailed scenario depends on the value 
of the response time r. Concerning the probability density in equilibrium, the attractor, and the conductivity in 
nonequilibrium, the properties of the standard Nose-Hoover thermostat in the periodic Lorentz gas are closer to the 
properties of the Gaussian thermostat than to the properties of thermostating by deterministic scattering, although 
the Nose-Hoover thermostat and thermostating by deterministic scattering share the property of keeping the energy of 
the particle on average constant in nonequilibrium, and for both thermostats the probability densities for the velocity 
components are related to a canonical probability density. 

Concerning the bifurcation diagrams we find that the dynamics of the three variations of the Nose-Hoover ther- 
mostat are in general "more chaotic". Even for higher field strengths there exist pronounced chaotic regions. The 
separate control of p\ and p y leads to a phase space filling attractor up to high field strengths in nonequilibrium. The 
equilibrium properties and the bifurcation diagram in nonequilibrium of this model are qualitatively closest to the 
properties of thermostating by deterministic scattering in comparison to the other versions. 

For the Nose-Hoover thermostat with field dependent coupling to the reservoir the thermodynamic entropy pro- 
duction is gcncrically not equal to the phase space volume contraction in contrast to the other models. 

An important question would be to look for common properties of all deterministic and time reversible thermostats. 
So far, only the existence of a fractal attractor in nonequilibrium appears to be typical. A more detailed investigation 
of these thermostating mechanisms should in particular involve a quantitative comparison, as, e. g., by means of 
computing Lyapunov exponents. An answer to this question would be helpful for obtaining a general characterization 
of nonequilibrium steady states. 
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FIG. 1. (a) Elementary cell of the periodic Lorentz gas on a triangular lattice, (b) Definition of the relevant variables. 
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FIG. 3. Probability density g(p x ) for s = 0.5: (a) Nose-Hoover thermostat with r 2 = 0.01 (dotted curve), r 2 = 1 (dashed 
curve), t 2 = 1000 (long dashed curve) and thermostating by deterministic scattering (solid curve), (b) Nose-Hoover thermostat 
with field dependent coupling to the reservoir with r 2 = 1 (dotted curve), additional control of < p 4 > with r 2 = 1 (dashed 
curve), control of < p 2 > and < p 2 > separately with r 2 = 0.1 and r 2 = 1000 (long dashed curve) and thermostating by 
deterministic scattering (solid curve). 



9 





0.28 

0.1 0.5 0.9 1.3 

£ 

FIG. 4. Conductivity a for the Nose-Hoover thermostat, (a) r 2 = 0.01, (b) r 2 = 1000. The numerical uncertainty of each 
point is less then symbol size. 
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FIG. 5. Poincare section of (/?, sin 7), as denned in Fig. |l] at the moment of the collision, for field strength e — 1: Nose-Hoover 
thermostat (a) r 2 = 0.01, (b) t 2 = 1000, (c) Nose-Hoover thermostat with field dependent coupling to the reservoir, r 2 = 1 
(d) additional control of p 4 , 
deterministic scattering. 



t = 1, (e) control of p x and p y separately, r x = 0.1 and r y = 1000, and (f) thermostating by 



11 




12 





2 4 6 8 10 



(b) 3 
2 

CO. 




























2 4 6 8 10 

8 

FIG. 7. bifurcation diagram for the dissipative limit (a) £ c = 1.0, (b) £ c = 2.0 
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FIG. 8. bifurcation diagram for the variation with field dependent coupling to the reservoir, r 2 = 1 
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FIG. 9. bifurcation diagram for the additional control of p 4 , r 2 = 1 
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FIG. 10. bifurcation diagram for the control of p\ and p 2 y separately, (a) r 2 = r 2 = 1, (b) rj = 0.1, r 2 = 1000 
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TABLE I. Phase space volume contraction rate and thermodynamic entropy production for the Nose-Hoover thermostat 
with field dependent coupling to the reservoir. The numerical error is < 0.001 



e 


r 2 = 1 




r 2 = 1000 






- < divf > 


Std 


- < divf > 


Std 


0.5 


0.152 


0.145 


0.145 


0.147 


1.0 


0.547 


0.561 


0.567 


0.592 


1.5 


1.240 


1.366 


1.256 


1.391 
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